Modulational Instability in Bose-Einstein Condensates under Feshbach Resonance Management

We investigate the modulational instability of nonlinear Schr{\"o}dinger equations with periodic variation of their coefficients. In particular, we focus on the case of the recently proposed, experimentally realizable protocol of Feshbach Resonance Management for Bose-Einstein condensates. We derive the corresponding linear stability equation analytically and we show that it can be reduced to a Kronig-Penney model, which allows the determination of the windows of instability. The results are tested numerically in the absence, as well as in the presence of the magnetic trapping potential

On the Modulational Instability of the Nonlinear Schr\"odinger Equation with Dissipation

The modulational instability of spatially uniform states in the nonlinear Schr\"odinger equation is examined in the presence of higher-order dissipation. The study is motivated by results on the effects of three-body recombination in Bose-Einstein condensates, as well as by the important recent work of Segur et al. on the effects of linear damping in NLS settings. We show how the presence of even the weakest possible dissipation suppresses the instability on a longer time scale. However, on a shorter scale, the instability growth may take place, and a corresponding generalization of the MI criterion is developed. The analytical results are corroborated by numerical simulations. The method is valid for any power-law dissipation form, including the constant dissipation as a special case

Soliton Dynamics in Linearly Coupled Discrete Nonlinear Schr\"odinger Equations

We study soliton dynamics in a system of two linearly coupled discrete nonlinear Schr\"odinger equations, which describe the dynamics of a two-component Bose gas, coupled by an electromagnetic field, and confined in a strong optical lattice. When the nonlinear coupling strengths are equal, we use a unitary transformation to remove the linear coupling terms, and show that the existing soliton solutions oscillate from one species to the other. When the nonlinear coupling strengths are different, the soliton dynamics is numerically investigated and the findings are compared to the results of an effective two-mode model. The case of two linearly coupled Ablowitz-Ladik equations is also investigated.Comment: to be published in Mathematics and Computers in Simulation, proceedings of the fifth IMACS International Conference on Nonlinear Evolution Equations and Wave Phenomena: Computation and Theory (Athens, Georgia - April 2007

Reaction-diffusion spatial modeling of COVID-19: Greece and Andalusia as case examples

We examine the spatial modeling of the outbreak of COVID-19 in two regions: the autonomous community of Andalusia in Spain and the mainland of Greece. We start with a 0D compartmental epidemiological model consisting of Susceptible, Exposed, Asymptomatic, (symptomatically) Infected, Hospitalized, Recovered, and deceased populations. We emphasize the importance of the viral latent period and the key role of an asymptomatic population. We optimize model parameters for both regions by comparing predictions to the cumulative number of infected and total number of deaths via minimizing the $\ell^2$ norm of the difference between predictions and observed data. We consider the sensitivity of model predictions on reasonable variations of model parameters and initial conditions, addressing issues of parameter identifiability. We model both pre-quarantine and post-quarantine evolution of the epidemic by a time-dependent change of the viral transmission rates that arises in response to containment measures. Subsequently, a spatially distributed version of the 0D model in the form of reaction-diffusion equations is developed. We consider that, after an initial localized seeding of the infection, its spread is governed by the diffusion (and 0D model "reactions") of the asymptomatic and symptomatically infected populations, which decrease with the imposed restrictive measures. We inserted the maps of the two regions, and we imported population-density data into COMSOL, which was subsequently used to solve numerically the model PDEs. Upon discussing how to adapt the 0D model to this spatial setting, we show that these models bear significant potential towards capturing both the well-mixed, 0D description and the spatial expansion of the pandemic in the two regions. Veins of potential refinement of the model assumptions towards future work are also explored.Comment: 28 pages, 16 figures and 2 movie

Healing Length and Bubble Formation in DNA

We have recently suggested that the probability for the formation of thermally activated DNA bubbles is, to a very good approximation, proportional to the number of soft AT pairs over a length L(n) that depend on the size $n$ of the bubble and on the temperature of the DNA. Here we clarify the physical interpretation of this length by relating it to the (healing) length that is required for the effect of a base-pair defect to become neligible. This provides a simple criteria to calculate L(n) for bubbles of arbitrary size and for any temperature of the DNA. We verify our findings by exact calculations of the equilibrium statistical properties of the Peyrard-Bishop-Dauxois model. Our method permits calculations of equilibrium thermal openings with several order of magnitude less numerical expense as compared with direct evaluations

Lengthscales and Cooperativity in DNA Bubble Formation

It appears that thermally activated DNA bubbles of different sizes play central roles in important genetic processes. Here we show that the probability for the formation of such bubbles is regulated by the number of soft AT pairs in specific regions with lengths which at physiological temperatures are of the order of (but not equal to) the size of the bubble. The analysis is based on the Peyrard- Bishop-Dauxois model, whose equilibrium statistical properties have been accurately calculated here with a transfer integral approach

Modulational and Parametric Instabilities of the Discrete Nonlinear Schr\"odinger Equation

We examine the modulational and parametric instabilities arising in a non-autonomous, discrete nonlinear Schr{\"o}dinger equation setting. The principal motivation for our study stems from the dynamics of Bose-Einstein condensates trapped in a deep optical lattice. We find that under periodic variations of the heights of the interwell barriers (or equivalently of the scattering length), additionally to the modulational instability, a window of parametric instability becomes available to the system. We explore this instability through multiple-scale analysis and identify it numerically. Its principal dynamical characteristic is that, typically, it develops over much larger times than the modulational instability, a feature that is qualitatively justified by comparison of the corresponding instability growth rates

Pattern Forming Dynamical Instabilities of Bose-Einstein Condensates: A Short Review

In this short topical review, we revisit a number of works on the pattern-forming dynamical instabilities of Bose-Einstein condensates in one- and two-dimensional settings. In particular, we illustrate the trapping conditions that allow the reduction of the three-dimensional, mean field description of the condensates (through the Gross-Pitaevskii equation) to such lower dimensional settings, as well as to lattice settings. We then go on to study the modulational instability in one dimension and the snaking/transverse instability in two dimensions as typical examples of long-wavelength perturbations that can destabilize the condensates and lead to the formation of patterns of coherent structures in them. Trains of solitons in one-dimension and vortex arrays in two-dimensions are prototypical examples of the resulting nonlinear waveforms, upon which we briefly touch at the end of this review.Comment: 28 pages, 9 figures, publishe

The role of mobility in the dynamics of the COVID-19 epidemic in Andalusia

Metapopulation models have been a popular tool for the study of epidemic spread over a network of highly populated nodes (cities, provinces, countries) and have been extensively used in the context of the ongoing COVID-19 pandemic. In the present work, we revisit such a model, bearing a particular case example in mind, namely that of the region of Andalusia in Spain during the period of the summer-fall of 2020 (i.e., between the first and second pandemic waves). Our aim is to consider the possibility of incorporation of mobility across the province nodes focusing on mobile-phone time dependent data, but also discussing the comparison for our case example with a gravity model, as well as with the dynamics in the absence of mobility. Our main finding is that mobility is key towards a quantitative understanding of the emergence of the second wave of the pandemic and that the most accurate way to capture it involves dynamic (rather than static) inclusion of time-dependent mobility matrices based on cell-phone data. Alternatives bearing no mobility are unable to capture the trends revealed by the data in the context of the metapopulation model considered herein